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Abstract 

We review recent work on understanding the next to leading order corrections to the classical 
fields that dominate the initial stages of a heavy ion collision. We have recently shown that the 
leading In 1/x divergences of these corrections to gluon multiplicities can be factorized into the 
JIMWLK evolution of the color charge density distributions. 



1. Introduction: Glass and Glasma 

At large energies (small x) the hadron or nucleus wavcfunction is characterized by a 
saturation scale Qs 3> A^^^^ arising from the strong nonlinear interactions of the color 
field. In the Color Glass Condensate (CGC) (for reviews see [112] ) framework the small x 
part of the hadron wavcfunction is described in terms of a classical Weizsacker- Williams 
(WW) field radiated by the hard, large x, sources. The color sources p are stochastic 
variables fluctuating according to a probability distribution I4^j^[/9(a;j^)], where y is the 
rapidity scale separating fast and slow partons [3. . 

The matter during the first fraction of a fermi in a collision of two such objects is 
what we refer to as the Glasma [A^ . The glasma configuration after the collision, at times 
< r < l/Qs, consists of longitudinal chromomagnetic and -electric field which depend 
on the transverse coordinate on a typical scale ~ 1/Qs- As the system expands the fields 
are diluted and can be treated as particles, forming the leading order (LO) production 
is the contribution that is computed in the numerically solving the classical Yang-Mills 
equations [5|6|7j . 

In the following we are concerned with the next to leading order (NLO) in g, h or, 
equivalently, loop corrections to this classical field. At NLO one can produce pairs of 
quarks (see Refs. |8I9I10) ) or gluons (real corrections) and one must take into account 
one loop corrections to the classical field (virtual corrections). We shall argue that these 
corrections have logarithmically divergent contributions, which must then be resummed 
into the renormalization group evolution of the sources M^[/9(a;^)] f ll|12j . 
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2. Factorization theorem 



It is perhaps useful to look first at the weak field limit of the CGC, where particle 
production can be computed using fcr-factorization ([13], see e.g. [M] for an application 
to heavy ion collisions). The leading order multiplicity is 
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d^pj^dy as pj 
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For the real part of the Leading Log correction to this result one must take the corre- 
sponding expression for double inclusive gluon production 

d^fc. 



dA^ 



1 1 



and integrate it over the phase space of the second gluon (q^, yq). Note that at leading 
log accuracy we have here taken the multi-Regge kinematical limit, assuming that the 
two produced gluons are far apart in rapidity (see e.g. [151). The integral over yq diverges 
linearly (this is the general behavior of the gg — > gg scattering amplitude in the high 
energy limit t fixed, s ~ —u — > oo). This divergence is compensated (to the appropriate 
order in a^) by the real part of the BFKL evolution equation for ify{k±i). 

In the fully nonlinear case of AA collisions the fcT-factorization 
is broken (see e.g. [5' 16]), and one must solve the equations 
of motion to all orders in the strong classical field. The 
analogue of the unintegrated parton distributions (py{k±) 
is the color charge density distribution M^j^[jo]. These are 
similar in the sense that they are not (complex) wavefunc- 
tions but (at least loosely speaking) real probability dis- 
tributions. Factorization can be understood as a statement 
that one has found a convenient set of degrees of freedom 
in which one can compute physical observable from only 
the diagonal elements of the density matrix of the incoming 
nuclei. The difference is that when in the dilute case these 
degrees of freedom are numbers of gluons with a given mo- 
mentum, in the nonlinear case the appropriate variable is 
the color charge density and the relevant evolution equation 

is JIMWLK, not BFKL. The kinematical situation, how- i:.- , tj j r . , 

' ' i'lg. 1. Production oi two gm- 

ever, remains the same. To produce a gluon at a very large ons: the integral over Ay is di- 
rapidity (or a contribution in the loop integral of the virtual vergent. 
contribution with a large k^) one must get a large -f -momentum from the right- moving 
source. Thus one is probing the source at a large fc^, i.e. small distances in x~ , and the 
result must involve Wy [p] at a larger rapidity. The underlying physical interpretation of 
factorization is that this fluctuation with a large fc+ requires such a long interval in x'^ to 
radiated that it must be produced well before and independently of the interaction with 
the other (left moving and thus localized in x+) source. The concrete task is then to show 
that when one computes the NLO corrections to a given observable in the Glasma, all 
the leading logarithmic divergences can be absorbed into the RG evolution of the sources 
with the same Hamiltonian that was derived by considering only the DIS process. This 
is the proof |11|12|17] of factorization that we will briefly describe in the following. 
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3. Deriving JIMWLK factorization 



Consider the single inclusive gluon multiplicity which is a sum of probabilities to 
produce n + 1 particles, with the phase space of the additional n must be integrated out 



dN ^ l_ 

^ n=0 



(3) 



Because we have a theory with external color sources of order p ^ 1/g, all insertions of 
source appear at the same order |18j . A calculation using the Schwinger-Keldysh formal- 
ism leads to the following results: At LO, the multiplicity is obtained from the retarded 
solution of classical field equations (here (...) includes the appropriate normalization 
and projection to physical polarizations) 



d^p 



(4) 



The NLO contribution includes the one loop correction to the classical field and the 
component of the Schwinger-Keldysh (SK) propagator in the background field 



d'^p 



(...) 



(5) 



Because of the SK index structure (also /3 
satisfies an equation of motion with a retarded 
boundary condition), one can express the prop- 
agation of a small fluctuation a'^ {x) above the 
past light cone S as a functional derivative 
T„ of the LO classical field A'^{x) with re- 
spect to its initial condition on S: a^{x) — 
J_^^^a{u) ■ TuA'^{x). This leads after some 
rearrangements to the expression for the NLO 
contribution to the multiplicity as a functional 
derivative operator acting on the leading order 
result: 




Fig. 2. The 1-loop one and two point functions 
in the background field. 
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Here the two point function below the light cone Q'^'^{u,v) = J ■ 



is bilinear in the small fluctuation field [x) satisfying the small fluctuation equation of 



'{k)e 



motion with and initial condition given by a plane wave limj.o^_ 
see Fig. [2] for a pictorial representation of this structure. 

The leading logarithmic contribution comes from the longitudinal component of the 
integral over fc, the momentum of the initial plane wave perturbation (and the corre- 
sponding momentum in the one loop source term for the equation of motion satisfied by 
/?). This LLog part of the functional derivative ^ operator turns out to be precisely 
equivalent to the sum of the JIMWLK Hamiltonians describing the RG evolution of 
the source distributions Wj,[/o]. The fact that no other terms with the same logarithmic 
divergences appear is the crucial result for factorization. The JIMWLK Hamiltonian 
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(7) 



is most naturally expressed in terms of Lie derivatives Da{x±) operating on the Wilson 
lines formed from the classical field (U{xj_) ~ Pexpj— ig / da;^ (l/V'^) p{xj^,x^)^ for 
the nucleus moving in the +z direction) in terms of which the kernel in Eq. ([T]) is 

nx^,y^) = ^ f ^-"^ - -/^^ - [t/..C/t C; c,t (8) 

Let us conclude by summarizing some important aspects of the JIMWLK factoriza- 
tion theorem of Refs. |11I12| . We are interested in the high energy kinematical limit, 
with transverse momenta Qs and the energy s large; the dilute limit is that of BFKL 
physics. The relevant degrees of freedom in this framework are the color charge densities 
p of the fast partons and the classical fields of the small x ones. The color charges p ^ 1/g 
are parametrically large, and thus the problem is inherently nonperturbative, but there 
is, however, a consistent weak coupling or loop expansion. The primary observables of 
interest are inclusive single and multigluon multiplicities (not cross sections for produc- 
ing a fixed number of particles), which leads to results that can be expressed in terms of 
retarded (and advanced) propagators. We express these propagators as functional deriva- 
tives with respect to the field on an initial surface on the light cone, which can then be 
mapped to the functional derivatives of the color fields of the individual nuclei appearing 
in the JIMWLK Hamihonian. 
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